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ABSTRACT 

The  prevailing  view  that  submesoscale  fluctuations  (horizontal  wavelengths  less  than  a  few  kilometers  and 
vertical  wavelengths  less  than  a  few  hundred  meters)  are  dominated  by  internal  gravity  waves  is  tested  by 
measuring  Enel's  potential  vorticity.  11  --  (/  +  V  x  \)-VB.  where  the  buoyancy  B  =  -ghp/p.,.  Unlike 
geostrophic  or  nonlinear  Ertei  vorticity -carry  ing  motions,  internal  waves  have  no  Ertei  vorticity  fluctuations. 
Velocity  and  temperature  profile  surveys  beside  Ampere  Seamount  reveal  appreciable  Ertei  enstrophy.  and  thus 
a  significant  non-internal-wave  component,  on  horizontal  wavelengths  of  6-15  km  and  vertical  wavelengths 
of  50-380  m.  The  twisting  terms  are  negligible  and  the  relative  vorticities  less  than  0  2/.  so  the  anomalies  are 
in  geostrophic  balance. 

It  is  unlikely  that  the  anomalies  arise  from  stirring  of  the  large-scale  isopycnal  gradients  of  stretching  and 
planetary  Ertei  vorticity  as  this  would  require  stirring  lengths  of  thousands  of  kilometers.  The  most  likely  source 
appears  to  be  forcing  at  the  seamount,  but  generation  by  ( i )  dissipative  3D  turbulence  in  the  pycnocline  or  ( ii ) 
detrainmem  of  the  winter  mixed  layer  cannot  be  absolutely  ruled  out.  It  remains  to  determine  whether  the 
coexistence  of  internal  wave  and  Ertei  vorticity-carrying  fluctuations  characterizes  smaller  scales  (X  s  50  m. 
\,i  s;  5  km)  in  the  deep  ocean  away  from  topography  as  well. 


1.  Introduction 

Ocean  fluctuations  with  vertical  wavelengths  less 
than  a  few  hundred  meters  and  horizontal  wavelengths 
less  than  a  few  kilometers  are  usually  interpreted  as 
internal  gravity  waves  (e.g..  Garrett  and  Munk  1979 ). 
However,  Holloway  ( 1983)  and  Muller  ( 1984)  have 
questioned  this  view^  pointing  out  that  geostrophic  and 
nonlinear  Ertei  vorticity-carrying  motions  can  coexist 
on  these  same  scales.  The  generalization  of  geostrophic 
motion  to  include  high  Rossby  number  |  f  |//and  high 
Froude  number  i  1/-I/7V  describes  the  part  of  the  flow 
that  carries  Ertei  vorticity  (Muller  1984).  Manifesta¬ 
tions  of  Ertei  vorticity-carrying  dynamics  have  been 
referred  to  as  blobs,  lenses,  blini,  pancake  eddies,  geo¬ 
strophic  turbulence,  and  stratified  turbulence  in  the 
literature.  These  motions  are  distinct  from  three-di¬ 
mensional  dissipative  turbulence  in  that  their  diapycnal 
motions  are  weak  and  they  do  not  mix  fluid  across 
isopycnals.  The  relative  contributions  of  internal  waves 
and  Ertei  vorticity-carrying  fluctuations  in  the  ocean 
may  have  important  consequences  for  understanding 
the  high-horizontal-wavenumber  end  of  the  potential 
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enstrophy  cascade,  horizontal  mixing,  and  the  pro¬ 
duction  of  dissipative  turbulence. 

The  challenge  for  observationalists  is  to  distinguish 
between  internal  wave  and  Ertei  vorticity-carrying 
fluctuations  in  the  sea.  These  two  kinds  of  motion  have 
several  distinguishing  properties.  Internal  waves  have 
intrinsic  frequencies  /<  oo  <  /V.  where  /is  the  Coriolis 
frequency  and  /V  the  buoyancy  frequency,  while  Ertei 
vorticity-carrying  fluctuations  are  subinertial  ( «  </). 
This  allows  identification  of  large-scale  fluctuations 
( which  are  little  affected  by  advective  nonlinearity )  us¬ 
ing  moored  time  series  measurements.  However, 
smaller  scales  become  increasingly  Doppler  smeared. 
For  example.  Briscoe  ( 1977).  Muller  etal.  ( 1 978),  and 
Eriksen  ( 1978)  could  not  identify  the  dynamics  of  fine- 
structure  in  the  internal  wave  continuum  band  (/<?  oi/.; 

<  N)  of  amplitude  2  cm  s  1  and  vertical  wavelength 

<  10  m.  This  finestructure  contained  more  kinetic  en¬ 
ergy  than  could  be  explained  by  the  observed  potential 
energy  assuming  linear  internal  wave  dynamics.  Kunze 
et  al.  ( 1990)  confirmed  this  result  and  reported  excess 
potential  energy  in  the  near-inertial  band  («/.  =  /). 
They  argued  that  this  could  be  explained  by  near-in¬ 
ertial  waves  being  Doppler  smeared  across  the  contin¬ 
uum  band  and  continuum  internal  waves  being  Dopp¬ 
ler  smeared  into  the  near-inertial  band  without  having 
to  invoke  non-internal  wave  dynamics.  This  conjecture 
is  supported  by  the  work  of  Sherman  and  Pinkel  ( 1991 ) 
who  transformed  CTD  profile  time  series  of  finescale 
strain  £-  (where  £  is  the  isopycnal  displacement)  into 
isopycnal-following  coordinates,  thereby  removing  the 
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effects  of  vertical  Doppler  shifting  and  aliasing.  In  this 
“semi-Lagrangian"  frame,  most  of  the  hnescale  strain 
fluctuations  were  attributable  to  low-frequency  internal 
waves.  Anderson  and  Pinkel  (  1493)  reach  the  same 
conclusion  for  hnescale  shear. 

As  its  name  suggests,  an  unambiguous  property  of 
Ertel  vorticity -carry  ing  flow  is  Fuel's  (  1942 )  potential 
vorticity,  which  will  be  defined  here  as  the  dot  product 
of  the  absolute  vorticity  ( 2L>  +  V  x  V )  and  the  gradient 
of  the  buoyancy  V(  B  +  />) 

II  =  (212  +  V  X  V).  V(B  +  />).  (  I  ) 

where  the  planetary  vorticity.  212  =  ( 0.  212  cos(  lat ).  212 
sin(  lat ))  =  (0.  /  cot(  lat )./').  and  the  buoyancy  B  +  h 
=  -,i !?>/>/ p„  has  been  split  into  a  large-scale  background 
B( : ).  which  depends  only  on  depth,  and  submesoscale 
anomalies  />(  v,  \ ,  r,  /).  Neglecting  both  the  horizontal 
gradients  of  vertical  velocity,  n\  and  it,,  and 
the  meridional  component  of  the  planetary  vorticity 
212  cos(  lat )  as  compared  to  the  vertical  shears.  uz  and 
t’_-.  (  I  )  can  be  expanded  into  a  background  plus  linear 
and  nonlinear  perturbations 

II  =  J\\-  +  //)_.  +  (i\  -  ur) A'2 

i - 1 

(i)  (ii)  (iii) 

+  ( t\  uv)h.  -  bxv-  +  bYii -.  (2) 

1 .  .  i 

(iv)  (v)  (vi) 


where  subscripts  denote  derivatives  and  A’2  =  Bz  is  the 
background  buoyancy  frequency  profile.  The  back¬ 
ground  Ertel  vorticity  (i)  involves  only  the  large-scale 
stratification  B(z)  and  the  planetary  vorticity  /. 
The  linear  perturbation  comprises  buoyancy-gradient 
anomalies  (vortex  stretching)  (ii)  and  relative  vorticity 
(iii).  Nonlinear  Ertel  vorticity  anomalies  arise  from 
coupling  of  vertical  buoyancy-gradient  anomalies  and 
relative  vorticity  (iv)  and  the  twisting  terms,  (v)  and 
(vi).  which  involve  horizontal  vorticities  (vertical 
shears)  and  horizontal  buoyancy  gradients. 

Term  (ii)  scales  as  the  strain  term  (iii)  as  the 
vorticity  Rossbv  number  R(  =  f//,  term  (iv)  as  R^:, 
while  terms  (v)  and  (vi)  scale  as  the  Froude  number 
squared  Fr2  =  l'2./ A'2  =  (it]  +  i-:)/A'2.  If  the  strain 
Rossbv  number  |f| //.  and  Froude  number  Vz/N 
arc  small,  the  motion  is  geostrophic.  Thermal  wind 
then  implies 


V  =  ft;; 


fb; 

A'2f 


_1_ 


PE 

KE’  *  ) 


where  ku  and  k-  are  horizontal  and  vertical  wavenum¬ 
bers.  and  PE  and  KE  are  the  available  potential  and 
horizontal  kinetic  energies  (Kunze  1993).  From  (3). 
it  is  clear  that  the  vortex-stretching  (ii)  and  relative 
vorticity  ( iii )  terms  have  the  same  sign  ( Hoskins  et  al. 
1985).  and  their  ratio  is  proportional  to  the  energy 
ratio  PE/KE  and  a  normalized  aspect  ratio  R/  =  (A 'k„j 


Jk;)2  (length-scale  Burger  number)  for  geostrophic 
flow. 

For  high  Rossbv  numbers  (A// >  |f//|  £>  I  ). 
earth's  rotation  is  unimportant  and  Ertel  vorticity- 
carry  ing  motion  can  be  described  as  cyclostrophic  or 
as  stratified  turbulence  (  Riley  et  al.  1981;  Lilly  1983). 
Muller  (  1984 )  coined  the  term  "vortical  mode"  to  de¬ 
scribe  Ertel  vorticitv-carry  ing  fluctuations  at  all  Rossbv 
and  Froude  numbers.  We  shall  refer  to  them  as  Enel 
vorticity-carry  ing  dynamics  in  this  paper. 

Internal  waves  have  no  Ertel  vorticity  fluctuations 
because  they  propagate  reversibly  and  the  Ertel  vorticity 
of  a  water  parcel  can  only  be  modified  by  irreversible 
processes  ( Pedlosky  1979;  Haynes  and  McIntyre  1987). 
For  example.  Lelong  and  Riley  (  1992)  have  shown 
that  weak  triad  wave-wave  and  wave-vortex  interac¬ 
tions  do  not  modify  Ertel  vorticity  .  Vortices  appear  in 
numerical  simulations  of  forced  "internal  waves"  (  Ri¬ 
ley  et  al.  1981;  Staquet  and  Riley  1989;  Herring  and 
Metais  1989)  only  if  the  forcing  directly  induces  Ertel 
vorticity  anomalies;  Herring  and  Metais  ( 1 989 )  found 
that  a  numerically  simulated  field  of  random  internal 
waves  without  Ertel  vorticity  forcing  did  not  develop 
Ertel  vorticitv-carrying  motions. 

Measuring  finescale  Ertel  vorticity  is  a  daunting  un¬ 
dertaking.  Relative  vorticity  and  nonlinear  twisting  can 
contribute  significantly  on  the  submesoscale  ( e.g..  Pol¬ 
lard  and  Regier  1992).  requiring  measurements  of 
horizontal  gradients  of  u,  t\  and  b,  which  are  difficult 
to  estimate  reliably  in  the  ocean.  Muller  et  al.  ( 1988) 
attempted  to  determine  the  linear  Ertel  vorticity  per¬ 
turbation.  /  A  2(£.  +  f//).  where  =  b-J X2,  using  the 
IWEX  trimooring  data  to  estimate  relative  vorticity  f 
and  horizor  4  divergence  V — from  which  they  deter¬ 
mined  the  strain. 


— for  horizontal  separations  of  5-1000  m.  Unfortu¬ 
nately,  a  three-point  array  cannot  distinguish  vorticity 
and  divergence  if  the  dominant  signal  is  from  scales  at 
or  smaller  than  the  array  as  is  the  case  in  the  ocean 
( Lien  and  Muller  1992b).  Furthermore,  for  horizontal 
separations  less  than  100  m.  internal  wave  isopycnal 
slopes  tilt  horizontal  vorticity  ( vertical  shear)  into  ver¬ 
tical  vorticity  enough  so  that  the  linear  Ertel  vorticity 
is  not  a  reasonable  approximation  for  the  total  Ertel 
vorticity.  Thus.  Muller  et  al.'s  ( 1988)  conclusion  that 
they  identified  Ertel  vorticity  finestructure  is  not  ten¬ 
able. 

In  this  paper,  we  will  look  for  submesoscale  Ertel 
vorticity  anomalies  in  two  velocity  and  temperature 
profile  surveys  collected  adjacent  to  Ampere  Seamount. 
Ampere  is  a  bank  rising  to  within  53  m  of  the  surface 
at  35°03.5'N.  I2°52.2’  W  in  the  eastern  North  Atlantic 
( Fig.  I  ).  Its  eastern  flank  is  relatively  steep  ( slope  0. 1 5; 
critical  reflection  frequency  7/)  and  western  flank  more 
gentle  (slope  0.03;  critical  frequency  1.3/).  This  site 


He;.  I.  Bathymetry  in  fathoms  for  the  eastern  North  Atlantic  showing  the  location  of  Ampere  Seamount  at  35°04'N.  I2°52'W.  The  summit  rises  to  within  53  m  of  the  surface. 
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Fig.  2.  Sampling  around  Ampere  Seamount  includes  an  XBT  box  survey  around  the  seamount, 
seven  drifter  deployments,  two  cross-shaped  XCP  surveys,  and  six  CTD  casts.  Bathymetry  is  in 
meters.  The  XCP  surveys  on  the  eastern  flank  contain  27  and  28  profiles  with  spacings  of  0.3-1 
km. 


was  chosen  because  it  was  thought  that  the  likelihood 
of  finding  an  Ertel  vorticity  signal  would  be  greater  in 
the  vicinity  of  topography  where  boundary  mixing  and 
flow  separation  might  produce  shed  eddies  in  the  wake 
(Armi  1978;  D’Asaro  1988b;  Geyer  and  Signell  1990; 
Signell  and  Geyer  1991 ).  We  also  examine  the  thermal 
wind  balance  and  the  ratio  of  horizontal  divergence  to 
relative  vorticity. 

Since  the  inception  of  this  project,  Lien  and  Muller 
( 1992a)  have  developed  a  complete  set  of  consistency 
relations  in  terms  of  vertical  and  horizontal  wavenum¬ 
bers  for  distinguishing  linear  internal  waves  and  geos- 
trophy.  A  companion  paper  (Kunze  1993)  investigates 
the  dependence  of  the  energy  ratio  ( or  energy  Burger 
number)  PE/KE  on  the  normalized  aspect  ratio  (or 
length-scale  Burger  number)  (N\:/f\/,)2  in  this  data. 
The  observations  are  compared  with  the  expected  re¬ 
lationships  for  internal  waves  and  geostrophy  (Kunze 
et  al.,  1990;  D’Asaro  and  Morehead  1991;  T  ien  and 
Muller  1992a).  The  theoretical  relations  are  used  to 


partition  energy  between  internal  waves  and  geostro¬ 
phy.  These  results  are  consistent  with  the  conclusions 
of  this  paper  in  that  geostrophic  and  internal  wave 
fluctuations  are  shown  to  coexist  on  the  same  scales 
next  to  Ampere  Seamount,  but  the  Burger  number  ap¬ 
proach  is  able  to  probe  smaller  horizontal  scales. 


2.  Measurements 

The  measurements  consist  of  an  XBT  box  survey 
60  km  on  a  side  centered  on  the  seamount,  seven  drifter 
deployments  within  10  km  of  the  summit,  two  cross¬ 
shaped  expendable  current  profilers  (XCP)  surveys, 
and  six  CTD  casts  on  the  eastern  flank  (Fig.  2).  The 
drifters  were  drogued  to  follow  the  average  flow  be¬ 
tween  100  and  200  m  depth  using  a  high-drag  line  in 
that  depth  range  and  a  low-drag  line  to  the  surface 
(Drever  1993,  personal  communication).  They  were 
tracked  during  6-10  September  1988  by  mooring  a 
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radar  transponder  on  the  summit,  then  steaming  the 
ship  from  one  drifter  to  the  next  to  get  range  and  bear¬ 
ing  of  the  transponder.  Range  is  accurate  to  ±200  m 
and  bearing  to  ±2°  (or  350  m  at  10-km  range).  By 
comparison,  Loran-C  positioning  at  Ampere  was  only 
good  to  one  kilometer. 

The  XBT  box  exhibited  no  large-scale  impinging 
geostrophic  flow  in  excess  of  the  expected  ±5  cm  s  1 
tidally  induced  noise.  However,  the  drifters  revealed 
subinertial  flow  in  the  near  field.  On  the  northern  flank 
of  the  seamount  and  extending  east  in  what  resembles 
an  interior  boundary  layer,  the  drifters  moved  eastward 
with  average  speeds  of  5-7  cm  s  1  ( Fig.  3 ).  Tidal  fluc¬ 
tuations  of  ±5  cm  s  1  were  also  primarily  in  the  east 
velocity  component.  On  the  eastern  flank,  motion  was 
weak  and  random  for  the  first  three  days,  suggesting  a 
stagnation  point.  Based  on  the  lateral  shear  of  >03  J 
on  the  northeastern  flank  of  Ampere  Seamount  (Fig. 
3),  flow  separation  and  eddy  shedding  were  thought 
most  likely  to  occur  there.  Consequently,  the  first  cross¬ 
shaped  finescale  XCP  survey  was  conducted  in  the 
near-field  shear  zone  on  the  afternoon  of  8  September. 
The  second  survey  was  made  ~  5  km  to  the  southeast 
on  the  afternoon  of  9  September  (Fig.  2).  Six  CTD 
casts  to  a  depth  of  2000  m  were  also  made  on  the  east¬ 
ern  flank  to  determine  the  local  T,  a„  relation.  During 
the  course  of  the  second  XCP  survey,  the  drifters  on 
the  eastern  flank  accelerated  to  the  south,  reaching 
speeds  of  ~  10  cm  s  '.  This  is  reflected  in  the  south¬ 
ward  drifter  velocities  on  the  southeast  flank  in  Fig.  3. 
The  cause  of  this  acceleration  is  unknown,  but  may  be 
related  to  seamount-trapped  waves  ( Brink  1990). 

The  XCP  surveys  contained  28  and  27  probes  with 
12-14  XCPs  in  each  leg  and  drop  spacings  of  0.3-1 
km.  The  sampling  sequence  for  both  surveys  was  from 
west  to  east,  then  from  south  to  north.  The  first  survey 
was  effectively  4  km  in  diameter  and  the  second  7  km. 
Each  survey  was  completed  in  under  three  hours  in  an 
effort  to  minimize  aliasing  of  the  spatial  gradients  by 
temporal  variability.  The  crossings  of  the  east-west  and 
north-south  legs  were  ~2  h  apart.  A  full  description 
of  the  XCP  sampling  and  data  can  be  found  in  Ken- 
nelly  et  al.  ( 1989 ). 

Expendable  current  profilers  measure  the  horizontal 
velocity  (u.  u)  electromagnetically  relative  to  an  un¬ 
known  but  depth-independent  constant  ( that  is,  they 
measure  the  baroclinic  flow )  and  temperature  with  an 
XBT  thermistor.  These  quantities  are  sampled  at  0.3-m 
intervals  from  the  surface  to  ~  1600-m  depth.  The  up¬ 
per  50  m  is  contaminated  by  temporally  aliased  surface 
swell  (Sanford  et  al.  1982)  and  will  therefore  be  ex¬ 
cluded  from  analysis.  In  the  pycnocline,  the  oceanic 
signal  typically  falls  below  instrument  noise  for  wave¬ 
lengths  smaller  than  10  m.  Therefore,  standard  pre¬ 
analysis  XCP  processing  smooths  the  raw  data  with  a 
6-m  triangular  window  every  3  m.  The  same  smoothing 
and  subsampling  was  applied  to  the  XBT  data.  The 
XCP  rms  velocity  error  is  ±0.4  cm  s  1  and  the  tem- 
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Fig.  3.  Drifter  velocities  over  100-200-m  depth  around  Ampere 
Bank  for  6-9  September.  On  the  northern  flank  and  extending  east, 
flow  is  eastward  at  5-7  cm  s'1.  On  the  eastern  flank,  flow  is  weak 
and  random  prior  to  9  September.  The  first,  or  near-field.  XCP  survey 
( Fig.  2 )  was  located  in  the  strong  shear  on  the  northeast  flank  between 
the  eastward  flow  and  the  “stagnation  point." 

perature  error  ±0.06°C.  As  we  will  be  confining  our 
attention  to  vertical  wavelengths  less  than  400  m,  the 
unknown  barotropic  velocity  is  of  no  concern. 

We  will  use  the  survey  data  to  compute  vertical  and 
horizontal  gradients.  The  velocity  error  will  result  in 
shear  errors  of  ±10  4  s  1  =  0.03N,,  over  50  m  in  the 
vertical  and  ±10  6  s  1  =  0.0 13/over  4  km  in  the  hor¬ 
izontal.  The  temperature  error  corresponds  to 
±10 '3oC  m  1  over  50  m  equivalent  to  10  h  s  2 
=  0.1  A'2,  in  A’2  and  ±10  5  °C  m'1  over  4  km  equiv¬ 
alent  to  10  *  s  2  =  0.05  J’N„  in  bx  given  the  observed 
db/dT  sc  -10~3  m  s  2  °C~'.  Thus,  as  compared  to  / 
=  8.3  x  10  5  s  1  and  the  150-600-m  depth-averaged 
buoyancy  frequency  N„  as  3.5  X  10'3  s'1,  instrument 
noise  in  velocity  contributes  smaller  errors  than  tem¬ 
perature  to  an  Ertel  vorticity  signal. 

3.  Estimating  Ertel  vorticity  II 

a.  Requirements 

Estimating  the  different  components  of  (2)  and 
identifying  Ertel  potential  vorticity  anomalies  requires 
determination  of  both  horizontal  and  vertical  gradients 
of  u.  v,  and  b.  This  can  be  accomplished  with  cross¬ 
shaped  XCP  surveys  provided  that  ( i )  there  is  a  tight 
T-S  relation  so  that  temperature  can  be  used  as  a  proxy 
for  buoyancy,  (ii)  a  valid  separation  can  be  made  be¬ 
tween  the  background  /V2  profile  and  the  perturbation 
buoyancy  gradient  b: ,  ( iii )  temporal  variability  during 
the  3  h  of  the  surveys  does  not  alias  the  horizontal 
gradients,  and  ( i  v )  the  surveys  are  horizontally  coherent 
so  that  horizontal  gradients  can  be  estimated. 
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Figure  4  displays  the  /  -  n„  relation  just  east  of  Am¬ 
pere  based  on  the  six  CTD  casts.  The  relation  is  tight 
for  temperatures  /'<  1()°C  (depths  greater  than  1250 
m)  and  for  1 1.7°  <  T  <  16.2°C  ( 150-600-m  depth), 
allowing  temperature  to  be  used  as  a  proxy  for  buoy¬ 
ancy  in  these  depth  ranges.  Temperature  is  not  a  re¬ 
liable  measure  of  buoyancy  in  the  salt-stratified  ther- 
mostad  between  depths  of  600  and  1200  m.  a  feature 
of  the  Mediterranean  salt  tongue,  or  in  the  upper  150 
m  (  T  >  16°C)  where  water  mass  perturbations  are 
evident.  Therefore,  the  temperature-buoyancy  relation 
can  be  reliably  estimated  only  for  150-600  and  1250- 
1500  m. 

The  background  buoyancy  frequency  profile  V2  =  B 
was  determined  from  the  mean  temperature  profile  of 
the  25  XBTs.  and  therefore  represents  a  60-km  scale 
average.  Anomalies  b:  are  relative  to  this  average  pro¬ 
file.  The  XBT  data  did  not  extend  below  750  m.  so  for 
the  remainder  of  this  paper,  we  shall  confine  our  at¬ 
tention  to  the  150-600  m  depth.  In  this  range.  />(  /  ) 
=  Tdb/dT.  with  dh/dT  =  -10  'ms  :  °C  1  based 
on  least-squares  fits  to  the  CTD  data  (Fig.  5).  The 
XBT-average  buoyancy  frequency  V(  r)  decreases  from 
3.6  X  10  's'  at  150-m  depth  to  3.3  X  10  's  1  at 
400  m.  then  increases  to  3.7  X  10  *  s  'at  600  m.  The 
average  buoyancy  frequency  over  150-600  m.  A'„ 
=  3.5  X  10  ’s',  will  be  used  to  normalize  the  Ertel 
vorticitv  terms. 
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Fl<i.  5.  The  least-squares  fit  relation  between  temperature  and 
buoyancy  used  to  compute  buoyancy  frequency  A ;  and  stretching 
6;  in  (2).  Bootstrap  standard  deviations  are  denoted  by  stippling. 
Where  the  stippling  extends  over  the  entire  range,  the  relation  is  not 
usable.  This  limits  buoyancy  estimates  to  depths  150-600  m  and 
1 250- 1 6(H)  m.  The  thick  dotted  curve  corresponds  to  a  least-squares 
cubic  lit  over  7  =  1 1 .7°—  1 6.20C'. 


Flo.  4.  Density  it,  versus  temperature  7  from  the  six  CTD  profiles 
collected  on  the  eastern  flank  of  Ampere.  The  7 -u,  relation  is  tight 
for  temperatures  less  than  I0.3°C  ( depths  greater  than  1 250  m)  and 
between  1 1 .7°  and  I6.2°C  ( 1 50-61  M)-m  depth ).  Temperature  cannot 
be  used  to  infer  buoyancy  in  the  stippled  ranges  between  650  and 
1200  m  and  above  150-m  depth. 


To  accurately  determine  horizontal  gradients,  tem¬ 
poral  variations  during  the  3-h  surveys  must  not  seri¬ 
ously  alias  the  horizontal  structure.  Figure  6  displays 
vertical  wavenumber  spectra  of  the  sums  and  differ¬ 
ences  of  velocity  profiles  taken  at  the  same  location 
(within  0.5  km)  but  2  h  apart.  If  temporal  variation 
was  small,  the  "difference"  spectrum  should  lie  well 
below  the  “sum"  spectrum.  For  wavelengths  larger  than 
50  m,  the  difference  spectrum  is  a  factor  of  3-4  smaller 
than  the  sum  spectrum,  commensurate  with  the  ex¬ 
pected  wavenumber-independent  change  for  a  Garrett- 
Munk  (GM)  internal  wave  field  (Munk  1981).  By 
comparison,  for  profile  pairs  taken  5  min  and  0.5  km 
apart  ( not  shown ).  the  difference  spectrum  lies  below 
the  sum  for  vertical  wavelengths  as  small  as  30  m. 
Therefore,  the  equipartition  of  energy  between  differ¬ 
ence  and  sum  spectra  at  wavelengths  smaller  than  50 
m  in  Fig.  6  is  due  to  temporal  variability.  Subsequent 
a..  Jysis  will  involve  making  survey  fits  that  filter  out 
the  smaller  horizontal  wavelengths  ( higher  frequencies ) 
responsible  for  much  of  this  temporal  variability.  Sim¬ 
ulations  using  the  GM  model,  low-pass  filtered  to  re¬ 
move  the  high  horizontal  wavenumbers  suppressed  by 
these  fits,  indicate  that  temporal  contamina.  on  of  the 
fits  will  be  a  factor  of  2-3  smaller  than  the  raw  differ¬ 
ences  in  Fig.  6.  Thus,  temporal  contamination  will  be 
about  1 0T  for  A-  >  50  m.  Vertical  wavelengths  smaller 
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A  r  <  0.5  km  At  =  0.08  days 


FlCi .  6.  Vertical  wavenumber  spectra  lor  kinetic  energy  in  the  sums  (  Z.  solid )  and  differences  ( A.  dotted )  of  profile  pairs  taken  at  the 
same  position  (within  0.5  km)  but  2  h  apart.  The  spectra  are  presented  in  both  log-log  (left)  and  variance-preserving  (right)  formats. 
Instrumental  noise  is  indicated  by  stippling  and  is  smaller  than  the  thickness  of  the  lines  except  at  high  wavenumbers  ( <  10  m)  in  the 
log-log  spectra.  For  wavelengths  A.  >  50  m.  the  sum  spectrum  exceeds  the  difference  by  factors  of  3-4.  At  smaller  wavelengths,  there  is 
equipartition  between  the  sum  and  difference  spectra. 


than  50  m  will  be  treated  as  unreliable  for  determining 
horizontal  gradients,  an  assumption  in  agreement  with 
the  noise  in  the  horizontal  gradient  estimates. 

To  illustrate  that  the  surveys  are  horizontally  co¬ 
herent.  Fig.  7  displays  the  east-west  section  of  east  ve¬ 
locity  //,  north  velocity  v.  and  vertical  displacement  £ 
=  (T  -  (T))/(T-)  from  survey  2.  Coherent  fluctu¬ 
ations  extend  across  all  or  most  of  the  section  down  to 
wavelengths  of  ~50  m. 

Because  the  surveys  are  horizontally  coherent,  we 
can  estimate  the  horizontal  gradients  d/dx  from  the 
east-west  legs  and  d/dy  from  the  north-south  legs  using 
linear  least-squares  fits  at  each  depth.  The  resulting 
gradients  can  be  combined  to  determine  dynamically 
relevant  variables.  By  fitting  over  many  points  in  the 
horizontal,  we  suppress  measurement  noise  and  the 
fluctuations  in  high  horizontal  wavenumber  that  con¬ 
taminated  the  IWEX  estimates  (Muller  et  al.  1988; 
Lien  and  Muller  1992b),  although  a  cross-shaped  array 
is  not  optimal  for  reducing  mutual  contamination  of 
relative  vorticity  and  horizontal  divergence  (Lien  1992. 
personal  communication;  see  appendix  B). 


Figure  8  displays  (a)  horizontal  divergence  (», 
+  t>,)//and  (b)  relative  vorticity  (vx  -  h,  )/  / profiles 
from  survey  1  for  several  horizontal  fitting  ranges  A r 
bracketing  survey  center.  As  the  fitting  scale  Ar  in¬ 
creases.  so  does  the  number  of  profiles  going  into  the 
fit  ( number  shown  at  base  of  each  profile ).  Therefore, 
the  standard  deviation  (stippling)  as  determined  by 
the  bootstrap  method  (Efron  and  Gong  1983)  dimin¬ 
ishes.  For  example,  at  Ar  =  4.3  km.  the  horizontal 
divergence  and  relative  vorticity  have  significant  fluc¬ 
tuations  down  to  wavelengths  of  50  m  with  rms  values 
of  0. 1-0.2/,  an  order  of  magnitude  larger  than  instru¬ 
ment  noise  but  only  a  factor  of  2  larger  than  unresolved 
oceanic  noise  ~0.07  /. 

Horizontal  buoyancy  gradients,  bx  and  bv,  were  es¬ 
timated  in  a  similar  manner  and  combined  with  vertical 
shears.  u:  and  v~.  averaged  over  fitting  scale  Ar,  to  ob¬ 
tain  the  twisting  terms  [(v)  and  (vi)  in  (2)].  This  ne¬ 
glects  contributions  from  subsurvey  scales  of  the  form 
((bK  -  (bx))( v:  -  <t';>)>. 

Inertially  backrotating  the  velocity  profiles  to  a 
common  time  before  carrying  out  the  fits  did  not  alter 
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FlCi.  7.  Sample  section  ol'(a)  east  velocity  u.  (b)  north  velocity  v,  and  (c)  vertical  displacement  (  - (T -  (  7'))/(  7.)  from  the  east- 
west  leg  of  the  second  XCP  survey.  Drop  numbers  are  indicated  along  the  bottom  axis.  The  dotted  vertical  lines  are  profile  zeros.  Structure 
down  to  vertical  wavelengths  of  ~50  m  can  (sc  traced  across  all  or  most  of  the  section. 


details  of  the  vorticity  and  divergence  profiles  appre¬ 
ciably  because  of  the  short  duration  of  the  surveys  ( 3 
h).  Likewise,  transforming  the  ».  v,  and  b  profiles  into 
isopycnal  coordinates  did  not  modify  the  estimates  be¬ 
cause  isopycnal  displacements  on  the  scale  of  the  sur¬ 
veys  were  small  compared  to  the  50-400-m  vertical 
wavelengths  of  interest. 

h  hr  tel  vorticity  profiles 

With  the  ability  to  estimate  relative  vorticity.  vertical 
shear,  and  buoyancy  gradients,  the  Ertel  vorticity  ( 2 ) 
can  be  assessed.  As  will  be  shown  (Fig.  10),  the  non¬ 
linear  terms  are  negligible.  Therefore,  profiles  of  the 
linear  perturbation  Ertel  vorticity.//)-  +  fiV2,  are  dis¬ 
played  in  Fig.  9.  The  vertical  buoyancy-gradient 
anomaly  b:  for  the  vortex-stretching  term  fbz  was  av¬ 
eraged  over  the  same  fitting  scale  as  relative  vorticity. 
The  linear  perturbation  Ertel  vorticity  has  significant 
anomalies.  For  example,  a  50-m  wavelength  feature 
with  amplitude  ~0.5  fN2  resides  at  400-m  depth  for 
Ar  =  4.3  km.  Smaller  amplitude  fluctuations  are  ap¬ 
parent  above  and  below  this  feature. 


c.  Spectra 

The  relative  contributions  of  the  various  linear  and 
nonlinear  terms  to  the  Ertel  vorticity  perturbation  ( 2 ) 
are  most  clearly  illustrated  with  vertical  wavenumber 
spectra.  Fourier  transforms  were  applied  to  de-meaned 
profiles  of  1 28  points  ( 380  m)  starting  at  1 50-m  depth 
and  windowed  at  both  ends  with  10%  sin2  tapers.  The 
spectra  are  displayed  in  Fig.  10.  Table  1  summarizes 
the  rms  contributions  to  the  Ertel  vorticity  for  vertical 
wavelengths  greater  than  50  m.  For  the  4-  and  7-km 
fits,  the  nonlinear  terms,  £b:,  bxv:,  and  bpc,  contribute 
little  and  are  below  oceanic  noise.  Since  the  Rossby 
number  is  less  than  0.2  and  the  nonlinear  terms  less 
than  0.05 //V,2,,  any  Ertel  vorticity  anomalies  can  be 
attributed  to  geostrophic  motion.  Likewise,  any  inter¬ 
nal  wave  fluctuations  will  be  linear. 

The  vortex-stretching//);  spectral  level  (thick  solid ) 
far  exceeds  that  of  the  other  Ertel  vorticity  terms.  In 
survey  1,  most  of  the  vortex-stretching  variance  (rms 
magnitude  0.2 fN2)  is  at  200-  and  50-m  wavelengths. 
Relative  vorticity  {N2  (thin  solid.  Fig.  10a)  also  has 
most  of  its  variance  ( rms  magnitude  0. 1  fN2)  at  these 
wavelengths  but  is  markedly  weaker. 
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Flo.  8.  Horizontal  divergence  V  =  (V- V„)  =  w,  +  vy  (a)  and  relative  vorticity  f  =  (V  X  V„)  =  vx  -  uy  (b)  as  functions  of  depth  and 
fitting  scale  Ar  based  on  horizontal  least-squares  linear  fits  to  the  velocity  profile  sections.  The  number  of  profiles  going  into  the  fit  is  marked 
along  the  bottom  axis  and  the  bootstrap  standard  deviations  about  the  fits  indicated  by  stippling.  At  fitting  scales  of  3. 5-4. 3  km.  horizontal 
divergence  and  relative  vorticity  have  significant  fluctuations  down  to  —  50-m  wavelengths  with  rms  values  of  0.2 /. 
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Fig  .  9.  Linear  perturbation  Ertel  vorticity.  fb.  +  fA2.  as  a  function  of  depth  and  fitting  scale  Ar  from  the  near-held  survey  ( left  panel ) 
in  the  upper  depth  range  where  the  T-a,  relation  is  tight.  The  background  Ertel  vorticity  /A2 //A’,2,  (right  panel)  is  based  on  the  mean  XBT 
profile.  The  dotted  vertical  lines  correspond  to  profile  zeros  and  are  located  at  the  appropriate  fitting  scale.  Bootstrap  standard  deviations 
are  indicated  by  stippling.  Estimates  are  not  plotted  if  the  standard  deviation  exceeds  1.5/A2. 
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Fir..  HI.  Vertical  wavenumber  spectra  of  the  two  linear  and  three  nonlinear  components  of  the  Ertel  vorticity  perturbation 
( normalized  by  IS:,.  ,V„  =  3.5  >  10  1  s  1 )  from  survey  I  for  a  fitting  scale  Ar  =  4.3  km  (a)  and  survey  2  for  a  fitting  scale 
Ar  -  6.8  km  (b).  Spectra  are  presented  in  both  log-log  (left  panels)  and  variance-preserving  (right  panels)  formats.  The 
shading  about  the  vortex-stretching  and  relative  vorticity  curves  denotes  one  standard  deviation.  Vortex  stretching / h.  ( thick 
solid)  dominates  except  at  the  lowest  resolved  vertical  wavenumber  (\=  -  380  m)  in  survey  2  where  vortex-stretching  and 
relative  vorticity  levels  are  comparable.  The  nonlinear  terms  are  negligible. 
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f  or  survey  2  ( fitting  scale  Ar  =  6.8  km  ).  the  vortex 
stretching  again  dominates  (Fig.  l()b).  Most  of  the 
vortex-stretching  variance  (rms  magnitude  0.25/.Y') 
is  at  380-.  I40-.  and  5()-m  wavelengths.  The  relative 
vorticity  is  much  weaker  at  I  10- and  50-m  wavelengths, 
indicating  an  Frtel  vorticity  perturbation  almost  en¬ 
tirely  associated  with  vortex  stretching.  However.  at 
780-m  wavelength,  the  relative  vorticity  and  vortc" 
stretching  variances  are  comparable,  so  the  signal  could 
be  due  to  internal  wave  or  geostrophic  motions. 

Which  of  these  is  responsible  can  be  determined  by 
examining  spectra  of  the  linear  Frtel  vorticity  anomaly  . 
/  />  t  f.Y-  (  Fig.  II).  For  linear  internal  waves,  vortex 
stretching  and  relative  vorticity  are  equal  and  opposite 
( /  b.  j'.Y: )  so  that,  if  the  relative  vorticity  was  largely 
from  internal  waves,  one  would  expect  the  variance  of 
/  />.  4-  s'.Y  to  be  reduced  from  that  of//)-.  In  contrast. 
/  b.  and  /Y:  are  in  phase  for  geostrophy  (  3 ).  so  geo¬ 
strophic  /  /).•  4  f.Y:  variance  will  exceed//).-  variance. 
Figure  I  lb  reveals  that  lb  +  </Y:  has  i>  renter  variance 
than  l  b:  at  A.  =  380  m  in  survey  2.  signify  ing  that  the 
relative  vorticity  at  this  wavelength  is  associated  with 
geostrophy.  F  rom  (3).  the  corresponding  horizontal 
wavelength  should  be  "15  km  or  rouehly  twice  the 
size  of  the  array.  ( Finding  horizontal  wavelengths  larger 
than  the  array  does  not  invalidate  the  analysis  since 
the  fits  include  contributions  from  all  wavelengths 
larger  than  the  array.)  The  vortex-stretching  peaks  at 
A  110  and  50  m  have  inferred  horizontal  wave¬ 
lengths  of  7  km  and  >4  km  where  the  presence  of 
some  internal  wave  vorticity  at  50-m  vertical  wave¬ 
length  leads  to  potential  underestimation  of  its  hori¬ 
zontal  wavelength.  Likewise,  in  survey  I  (  Fig.  I  la ).//>_- 
+  i'X:  exceeds  vortex  stretching  at  the  peaks,  signifying 
that  the  bulk  of  the  relative  vorticity  is  in  phase  with 
the  vortex  stretching  and  therefore  associated  with 
geostrophy.  The  corresponding  horizontal  wavelengths 
(  3 )  at  A..  =  200  and  60  m  are  1 2  and  >3  km.  respec¬ 
tively.  The  fractions  of  kinetic  energy  contributed  by 
internal  waves  and  geostrophy  are  considered  further 
in  section  3d  and  by  Kunze  (  1993). 


While  vortex  stretching  greatly  exceeds  CiM  model 
levels  (Fig.  1 2a ).  the  vorticity  spectra  are  roughly  com¬ 
parable  except  at  A.  ;  60-100  m  in  survey  1  and  A 
-  380  m  in  suivey  2  (  Fig.  12b).  For  this  comparison, 
the  high  hori/ontal  wavenumbers  suppressed  in  the 
measurements  by  the  least-squares  lilting  are  filtered 
from  the  CiM  model  with  a  simple  boxcar  low-pass 
filter  with  cutoff  at  k,,  =  7r/Ar(  appendix  A)(theGM 
spectra  differ  for  the  two  surveys  because  of  their  dif¬ 
ferent  cutoff  wavenumbers).  Given  the  well-known 
factor  of  three  variability  in  energy  density  on  these 
scales  (  Munk  1981  ).  more  elaborate  attempts  to  lit  the 
GM  model  to  the  data  arc  not  justified. 

d  Relative  vorticity  f  and  horizontal  divergence  Vt/ 

Since  the  twisting  terms  are  small,  so  that  isopycnal 
tilting  can  be  neglected,  relative  vorticity  j"  =  V  x  V/ 
and  horizontal  divergence  V  =  V  •  X,,  (  =  -n_  -  -  £., ) 
can  also  be  used  to  distinguish  linear  internal  waves 
and  geostrophy  since.  ( i )  for  geostrophic  motions,  hor¬ 
izontal  divergence  is  weak  compared  to  vorticity  (V 
<1  f).  while  for  internal  waves.  <r?>/(  f:>  =  u //  / : 
(Fofonolf  1969).  so  that  ( ii )  near-inertial  waves  have 
comparable  divergence  and  vorticity  (V  ^  (')  and  ( iii ) 
high-frequency  internal  waves  have  divergence  ex¬ 
ceeding  vorticity  (V  {').  Thus,  nonnegligible  diver¬ 
gence  indicates  the  presence  of  internal  waves,  and 
vorticity  in  excess  of  the  divergence  is  ascribable  to 
geostrophic  motions.  Relative  vorticity  anq  .orizontal 
divergence  are  computed  in  the  same  way.  This  reduces 
concerns  about  different  sampling  and  processing  that 
arise  when  combining  horizontal  gradients  of  it  and  r 
with  vertical  gradients  of  buoyancy. 

The  observed  relative  vorticity  and  horizontal  di¬ 
vergence  fields  are  compared  with  GM  model  spectra 
for  these  quantities  in  Fig.  13.  High  horizontal  wave- 
numbers.  suppressed  in  the  observations  by  fitting  ovu 
the  surveys,  have  been  removed  from  the  GM  spectra 
with  a  boxcar  filter  (appendix  A  )  for  this  comparison. 
Filtering  has  its  greatest  impact  on  the  horizontal  di¬ 
vergence.  which  is  dominated  by  high-frequency,  high- 
horizontal-wavenumber  fluctuations  in  the  GM  model. 
More  sophisticated  filtering  of  the  GM  model  to  ac¬ 
count  for  contamination  of  the  vorticity  by  the  diver¬ 
gence  (Lien  and  Miillcr  1992b)  is  not  warranted  be¬ 
cause  this  effect  boosts  the  GM  vorticity  level  by  at 
most  I  O'?  above  the  divergence  for  vertical  wavelengths 
A-  <  60  m  (Lien  1992.  personal  communication,  see 
appendix  B). 

For  the  near-field  survey  I  (Fig.  13a).  the  observa¬ 
tions  share  with  the  low-pass  horizontal-wavenumber- 
filtered  GM  model  the  tendency  for  ( i )  divergence  to 
exceed  vorticity  at  low  vertical  wavenumbers,  and  ( ii ) 
divergence  and  vorticity  to  be  comparable  at  high  ver¬ 
tical  wavenumbers.  However,  the  relative  vorticity 
variance  is  twice  the  horizontal  divergence  variance 
for  A;  =  60-100  m.  Taking  the  vorticity  in  excess  of 
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Fl(i.  12.  Comparison  of  observed  vertical  wavenumber  spectra  for  vortex  stretching  J'b.  (a)  and  relative  vorticity  f.V*  (b) 
from  survey  I  <  Ar  4.3  km:  thick  solid)  and  survey  2  ( Ar  =  6.8  km;  thin  solid)  with  the  GM  model  spectra  low-pass  filtered 
to  remove  the  high  horizontal  wavenumbers  (A;/  >  *7  Ar)  suppressed  by  fitting  over  the  surveys  (dotted).  All  spectra  have 
been  normalized  by  / ,Y;,.  The  observed  vortex-stretching  spectra  ( a )  greatly  exceed  GM.  while  the  observed  relative  vorticity 
spectra  ( b )  have  comparable  levels  and  shapes  to  the  GM  model  except  for  excess  vorticity  at  60- 1 00  m  wavelength  in  survey 
I  and  at  the  largest  resolved  wavelength  ( A.  =  380  m)  in  survey  2. 


JOl'RNAI  Ot  PHYSICAL  (H  I  ANOGR  APHY 


A  =  4.3  km 


100  90  30  10  n 


300  100  SO  30  10m  _ (a) 

— i - 1 - 1 - 1 - 1 - 1  0.08 


v/f 

GM  V/f 

C/f 

GMC/f 


k  (rad/m) 


100  SO  30 


10  1 
(rad/m) 


A  r2  =  6.8  km 


v/f 

GM  v/f 

C/f 

GM£/f 


10  1 

kz  (rad/m) 


k  (rad/m) 


Fig.  13.  Comparison  of  the  measured  vertical  wavenumber  spectra  for  horizontal  divergence  (V- V„)//(thick  solid)  and 
relative  vorticity  (V  x  V„)//(thin  solid)  from  survey  I  for  Ar  =  4.3  km  (a)  and  from  survey  2  for  Ar  --  6.8  km  (b).  with 
GM  model  spectra  low-pass  filtered  to  remove  high  horizontal  wavenumbers  (k„  >  jr/ Ar)  suppressed  by  fitting  over  the 
surveys  ( thick  and  thin  dotted ).  Bootstrap  error  bars  (stippling  (are  less  than  the  width  of  the  curves  in  the  variance-preserving 
format  ( right  panels).  Vorticity  exceeds  divergence  consistent  with  gcostrophic  but  not  internal  wave  dynamics  at  A-  =  50- 
1  (X)  m  in  survey  I  and  at  the  largest  resolved  wavelength  ( Ar  =  380  m)  in  survey  2. 
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the  divergence  to  he  geostrophic.  this  suggests  that  at 
least  one-third  of  the  kinetic  energy  in  this  hand  is  as¬ 
sociated  with  geostrophic  motion.  This  is  commen¬ 
surate  with  the  conclusion  from  Tig.  I  la  that  the  rel- 
ative  sorticity  in  this  wavenumber  hand  was  largely 
geostrophic.  Likewise,  for  survey  2  (  Fig.  1 3b).  relative 
vorticity  exceeds  horizontal  divergence  by  a  factor  of 
4  at  A-  380  m.  signify  ing  0.5-4  times  as  much  kinetic 
energy  in  geostrophic  as  internal  wave  motions  at  this 
wavelength.  I  his  is  in  agreement  with  the  conclusion 
from  Fig.  lib.  For  vertical  wavelengths  between  50 
and  200  m.  geostrophic  kinetic  energy  is  0.05-2  times 
that  in  internal  waves  because  the  vorticity  variance  is 
twice  the  divergence  variance. 

e.  'I'cstinfi  thermal  wind 

Ft  is  common  to  test  for  geostrophy  by  comparing 
the  two  terms  of  the  thermal  wind  relation. 

/  V- =  ( -t  -.  ».)  and  V ///>  =  (/),. />,  ). 

For  geostrophic  flow,  the  ratio  (./^L  :)/((V„/>)2) 
=  I.  while  for  internal  waves. 

(/'  F  r )  =  /:A;KE  _  /2(  A'2  -  ar):(w2  4-  / 2 ) 

<(V„/>):)  A2Aj/PE  A -4(w2-/2)2 

(\2kh+2f2k2;)f2k*; 

A'4A//A 2  4 

( Fofonoff  1969).  The  internal  wave  ratio  exceeds  one 
for  ce  <  1.7/and  is  less  than  unity  for  higher  frequen¬ 
cies.  Measured  vertical  wavenumber  spectra  for  vertical 
shear/ V.-,  horizontal  buoyancy  gradient  V,,/?.  and  the 
ageostrophic  acceleration  V.,  =  (  -  A,  -/>,  -fu:) 
are  compared  with  GM  spectra  for  these  quantities  in 
Fig.  14.  The  GM  spectra  have  been  low-pass  filtered 
in  horizontal  wavenumber  to  remove  the  high  hori¬ 
zontal  wavenumbers  suppressed  by  fitting  over  the  sur¬ 
veys,  and  high-pass  filtered  in  frequency  to  eliminate 
frequencies  below  1.004/  that  is,  below  the  near-in¬ 
ertial  peak  (D'Asaro  and  Perkins  1984). 

If  the  fluctuations  were  purely  geostrophic.  the  ob¬ 
served  /I and  V//  spectra  would  be  identical  and  the 
ageostrophic  acceleration  I  spectra  would  vanish. 
While  the  measured  horizontal  buoyancy-gradient 
spectra  contains  more  variance  than  the  GM  model, 
the  vertical  shear  contains  much  more  variance  than 
the  horizontal  buoyancy  gradient  for  both  surveys,  sig¬ 
nifying  a  large  near-inertial  wave  contribution. 

This  apparent  contradiction  with  the  conclusions  of 
Figs.  9-13  can  be  reconciled  by  recognizing  that  the 
velocity  and  buoyancy  fits  contain  contributions  from 
horizontal  scales  at  and  larger  than  the  surveys.  If  the 
horizontal  wavenumber  spectra  are  red  with  slopes 
more  positive  than  -3.  then  the  observed  vortex 
stretching  J'h.  and  vertical  shear/!  - will  be  dominated 
by  the  energy-containing  horizontal  scales,  which  could 


be  larger  than  the  survey,  while  the  dominant  contri¬ 
butions  to  relative  vorticity  f.V2  and  horizontal  buoy¬ 
ancy  gradient  V Hh  will  be  from  scales  comparable  to 
the  surv  ey  size.  This  highlights  the  limitations  of  using 
thermal  wind  to  test  for  geostrophy  in  surveys  when 
broadband  horizontal  internal  wave  fluctuations  co¬ 
exist  on  the  same  scales  as  the  geostrophic  flow.  This 
problem  does  not  exist  for  the  Ertel  vorticity  test  (  Figs. 
9-1  I ).  which  isolates  the  Ertel  vorticity -carry  ing  part 
of  the  flow  from  internal  w  av  es. 

4.  Summary 

Velocity  and  temperature  profile  surveys  near  Am¬ 
pere  Seamount  have  been  used  to  examine  three  pairs 
of  consistency  relations  for  internal  waves  and  geos¬ 
trophy  ( Table  2  ):  the  Ertel  potential  vorticity  .  the  ratio 
of  horizontal  divergence  to  relative  vorticity.  and  ther¬ 
mal  wind.  A  fourth  consistency  relation  pair  involving 
the  ratio  of  potential  to  kinetic  energy  is  investigated 
in  a  companion  paper  ( Kunze  1993).  Although  unable 
to  shed  light  on  the  controversial  dy  namics  of  fluctu¬ 
ations  with  vertical  wavelengths  less  than  50  m  (Muller 
et  al.  1988;  Kunze  et  al.  1990;  Sherman  and  Pinkel 
1991;  D'Asaro  and  Morehead  1991;  Lien  and  Muller 
1991;  Anderson  and  Pinkel  1993 ).  the  analysis  reveals 
submesoscale  Ertel  vorticity  anomalies  with  rms  mag¬ 
nitudes  of0.2-0.4/V2(  Table  I  )on  vertical  wavelengths 
as  small  as  50  m  ( Figs.  9  and  10).  scales  traditionally 
thought  to  be  dominated  by  internal  waves.  These 
anomalies  are  an  indisputable  signature  of  geostrophic 
or  nonlinear  Ertel  vorticity-carrying  dynamics  because 
internal  waves  and  seamount-trapped  waves  ( Brink 
1990;  Codiga  1993 )  have  no  Ertel  vorticity  fluctuations 
in  the  water  column.  Thus.  Ertel  vorticity-carrying 
dynamics  coexists  on  the  same  scales  as  internal  waves 
(at  least  in  the  vicinity  of  Ampere  Seamount).  Relative 
vorticities  were  less  than  0.2/and  the  nonlinear  Ertel 
vorticity  terms  were  weak  compared  to  the  linear  terms 
(Fig.  10;  Table  1 ):  therefore,  the  anomalies  are  in  geo¬ 
strophic  balance. 

As  on  basin  scales,  the  Ertel  vorticity  signal  over 
4-7  km  is  dominated  by  vortex  stretching  fh:  ( Fig. 
10):  the  relative  vorticity  /V2  was  much  weaker.  The 
vortex-stretching  signal  is  reminiscent  of  the  density 
layering  discussed  in  the  internal  wave/finestructure 
literature  (Phillips  1971;  Garrett  and  Munk  1971; 
McKean  1974:  Eriksen  1978:  Levine  and  Irish  1981 ). 
If  the  stretching  were  due  to  internal  waves,  it  would 
be  compensated  by  the  vorticity  (//>_-  =  -  /V2).  Since 
the  vortex  stretching  dominates,  it  can  be  concluded 
that  it  is  principally  due  to  geostrophy.  Moreover, 
where  the  vortex-stretching  signal  at  a  vertical  wave- 
number  is  attributable  to  geostrophy.  so  is  the  available 
potential  energy  since  PE  =  b]/(2k]N2).  This  is  com¬ 
mensurate  with  the  partition  of  energy  between  internal 
waves  and  geostrophy  found  using  the  energy  and 
length  scale  Burger  numbers  ( Kunze  1993).  From  this. 
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Fig.  14.  Comparison  of  vertical  wavenumber  spectra  of  vertical  shear/!',  (solid),  horizontal  buoyancy  gradient  V„/>  = 
(dotted),  and  the  ageostrophic  acceleration  V„  =  (  -bx  +  fv.,  -b,  -  fit, )  (dashed)  with  their  GM  model  spectra  low-pass 
filtered  to  remove  horizontal  wavenumbers  suppressed  by  the  fitting  and  high-pass  filtered  to  remove  frequencies  lower  than 
1 .004  /.  All  spectra  have  been  normalized  by  /  A’„.  If  the  observed  fluctuations  were  purely  geostrophic.  the  vertical  shear  and 
horizontal  buoyancy  gradient  would  be  identical  and  the  ageostrophic  acceleration  would  vanish.  The  shear  exceeds  the 
buoyancy  gradient,  and  the  ageostrophic  acceleration  is  of  comparable  magnitude  to  the  shear  in  the  observations,  qualitatively 
consistent  with  the  GM  model. 
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Tabi  i  2  Consistency  relations  for  internal  waxe  and  geostrophic 
fluctuations  examined  in  this  paper  including  the  ratio  of  vortex 
stretching  Ih.  to  relative  vorticily  f.V\  the  ratio  of  horizontal  diver¬ 
gence  C  to  relative  vorticity  and  the  ratio  of  vertical  shear  / 1  ;  to 
horizontal  buoyancy -gradient  V„/>  terms  in  the  thermal  wind  balance. 
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it  might  be  thought  that  vortex  stretching  alone  could 
be  used  to  estimate  Ertel  vorticity  down  to  scales  of  a 
few  kilometers.  However,  given  the  uniqueness  of  the 
measurement  site  and  that  internal  waves  also  con¬ 
tribute  to  submesoscale  vortex  stretching,  relative  vor¬ 
ticity  must  be  measured  to  ensure  that  fh:  is  not  com¬ 
pensated. 

Spectra  of  the  linear  Ertel  vorticity  anomaly,  fh: 
+  f  A’2  ( Fig.  11),  contain  more  variance  than  do  those 
of  vortex  stretching  or  relative  vorticity  alone.  This 
suggests  that  the  relative  vorticity  also  is  largely  due 
to  geostrophy  since  internal  wave-induced  vortex 
stretching  and  relative  vorticity  will  balance  each  other 
(  I fb:  +  fA’2  |  <  |  fb;  |  =  |  fAr2  J )  when  nonlinear  terms 
are  small.  Using  (3),  the  horizontal  wavelengths  for 
the  peak  bands  are  7-15  km.  The  geostrophic  inter¬ 
pretation  is  borne  out  by  comparison  of  vertical  wave- 
number  spectra  of  relative  vorticity  and  horizontal  di¬ 
vergence  (Fig.  13).  which  reveal  vorticity  in  excess  of 
divergence  in  peak  bands,  consistent  with  geostrophy 
but  not  internal  waves ( section  3d).  Divergence  dom¬ 
inates  in  other  bands,  consistent  with  internal  wave 
but  not  geostrophic  dynamics.  Thus,  there  are  contri¬ 
butions  to  the  kinetic  energy  from  both  internal  waves 
and  geostrophy. 

Using  velocity  profile  surveys  collected  in  the  Beau¬ 
fort  Sea.  D' Asaro  and  Morehead  (1991)  also  reported 
evidence  of  kinetic  energy  associated  with  Ertel  vortic- 
ity-carrying  dynamics  on  horizontal  scales  less  than  a 
few  kilometers  at  =  40  m.  At  0.2  cm2  s  2,  their 
signal  is  comparable  to  the  geostrophic  kinetic  energy 
found  beside  Ampere  (Kunze  1993)  and  three  times 
greater  than  that  of  high-frequency  internal  waves  in 
the  weak  Arctic  internal  wave  regime.  This  part  of  the 
Arctic  is  known  to  contain  many  small  eddies  ( Manley 
and  Hunkins  1985;  D’ Asaro  1988a).  which  D’ Asaro 
( 1988b)  argues  are  generated  by  frictional  torques  in 
Barrow  Canyon. 

In  contrast,  Kunze  et  al.  ( 1990),  Sherman  and  Pinkel 
( 1991 ),  and  Anderson  and  Pinkel  ( 1993)  concluded 
that  finescale  ( A- «  10  m)  shear  is  dominated  by  near- 
inertial  waves  based  on  measurements  1000  km  off  the 


coast  of  California.  Lien  and  Muller  (  1992b)  find  that 
fluctuations  with  A.-  >  70  m  are  consistent  with  internal 
waves  in  the  Sargasso  Sea.  Thus,  it  is  too  soon  to  tell 
whether  the  results  found  here  are  typical  of  the  ocean 
or  due  to  topography. 

5.  Discussion 

The  finding  of  strong  Ertel  vorticity  anomalies 
—0(0.3 /A’2)  beside  Ampere  Seamount  raises  the 
questions  of  how  they  were  generated  and  whether  they 
are  ty  pical  of  the  ocean  pycnocline  or  due  to  the  prox¬ 
imity  of  topography  (Fig.  1 ).  Four  formation  mecha¬ 
nisms  seem  plausible;  (i)  transfer  of  variance  from 
lower  horizontal  wavenumbers  as  part  of  the  potential 
enstrophv  cascade  to  high  k,,  of  geostrophic  turbulence 
(Charney  197 1 ),  ( ii )  dissipative  turbulence  and  mixing 
in  the  pycnocline  ( Riley  et  al.  1 98 1 ).  ( iii )  detrainment 
of  the  late-winter  mixed  layer  followed  by  subduction 
(e.g..  Stommel  1979;  Marshall  et  al.  1993),  and  (iv) 
intensified  mixing  processes  on  topography.  These  are 
discussed  below  in  order  of  increasing  plausibility. 

If  interpreted  as  geostrophic  turbulence,  the  Ertel 
vorticity  anomalies  must  be  produced  by  eddy  stirring 
of  the  large-scale  background  gradients  of /A,2_.  This 
implies  turbulent  stirring  lengths  5I1/|V///A'2|  of 
thousands  of  kilometers  based  on  the  rms  Ertel  vorticity 
anomaly  50  =  0.3 /TV2  and  the  large-scale  background 
gradients  of /A'2  due  to  the  /?  effect  and  gyre-scale  cir¬ 
culation  (McDowell  et  al.  1982;Keffer  1985 ).  The  ob¬ 
servations  occupy  a„  =  26.7-27.1  ( Fig.  4 ),  close  to  the 
homogenized  pool  described  by  McDowell  et  al. 
( 1982).  where  Rhines  and  Young  ( 1982)  predict  van¬ 
ishing  eddy  potential  energy;  however,  stronger  /TV2 
gradients  ~0(dA’,2,)  exist  to  the  north.  Stirring  lengths 
of  thousands  of  kilometers  are  unrealistically  large,  in¬ 
dicating  that  the  observed  signal  is  not  representative 
of  the  high-wavenumber  end  of  the  potential  enstrophv 
cascade. 

Because  turbulent  dissipation  and  mixing  are  weak 
and  confined  to  vertical  scales  of  a  few  meters  in  the 
ocean  interior,  Ertel  vorticity  anomalies  produced  by 
dissipative  turbulence  in  the  pycnocline  should  likewise 
be  weak  and  confined  to  small  vertical  scales.  Given 
typical  dissipation  rates  t  =  10  9  W  kg”1  (Mourn  and 
Osborn  1986;  Gregg  1987;  Gregg  1989;  Yamazaki  et 
al.  1990),  mixing  efficiencies  less  than  0.2  (Itsweire  et 
al.  1986),  and  assuming  equal  generation  of  internal 
waves  and  Ertel  vorticity-carrying  fluctuations,  an  up¬ 
per  bound  for  the  rate  of  energy  transfer  to  Ertel  vor¬ 
ticity  anomalies  is  ~0.3  cm2  s  2  day  "1.  This  is  suffi¬ 
cient  to  fill  the  —  1  cm2  s” 2  in  the  Ertel  vorticity  anom¬ 
alies  (Kunze  1993)  in  under  a  week  (lower  bound). 
However,  to  do  so  would  require  many  mixing  events, 
each  of  only  a  few  buoyancy  periods  duration,  to  be 
cumulative.  While  Gregg  et  al.  ( 1986)  found  that  one 
turbulent  patch  associated  with  a  near-inertial  wave 
packet  persisted  for  many  days,  most  turbulence  events 
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appear  to  occur  randomly  throughout  the  water 
column. 

Random  turbulence  will  erode  Ertel  vortieity  anom¬ 
alies  rather  than  create  them.  (Erosion  times  will  be 
much  longer  than  the  above-quoted  formation  times 
because  most  of  the  turbulent  dissipation  is  supplied 
by  small-scale  internal  waves.)  An  eddy  viscosity  of  r 
=  i/l  t  *£  t/(4A:)  =  0.1  X  10  4  m:  s  1  can  be  de¬ 
duced  assuming  that  internal  wave-induced  turbulence 
acts  to  smooth  larger-scale  shear  as  well  as  the  unstable 
small-scale  shear,  this  is  a  reasonable  assumption  if 
turbulence  production  is  driven  by  internal  wave-wave 
interactions  (Henyey  et  al.  1986)  uncorrelated  with 
the  Ertel  vortieity-carrying  flow.  The  above  eddy  vis¬ 
cosity  implies  a  decay  time  (Gill  1981  )  of  10  yr  for 

=  400  m  and  2-3  months  for  A.  =  50  m.  Corre¬ 
sponding  advection  length  scales  are  20  000  and  400 
km  for  a  5-cm  s  1  current.  Thus,  the  observed  50-m 
vertical  wavelength  feature  (Figs.  9  and  10)  will  not 
be  long-lived  or  of  distant  origin. 

Double  diffusion  can  also  create  stretching  vortieity 
in  the  form  of  staircases  and  thermohaline  intrusions 
but  the  associated  horizontal  velocity  holds  are  much 
weaker  ( <  1  cm  s  1 )  than  those  observed. 

An  ironclad  case  cannot  be  made  against  production 
by  microscale  dissipative  turbulence  in  the  pycnocline. 
However,  given  that  the  strongest  turbulent  dissipation 
and  mixing  occur  at  the  surface  (Schmitt  and  Olson 
1985:  Shay  and  Gregg  1986)  and  benthic  boundaries 
(Nabatovand  Ozmidov  1988:  Toole  and  Schmitt  1992. 
personal  communication),  these  might  be  argued  to 
be  the  most  likely  sites  for  generation  of  submesoscale 
Ertel  vortieity  perturbations. 

The  retreat  of  the  winter  mixed  layer  is  likely  to 
leave  behind  “granularity"  in  Ertel  vortieity  just  as  it 
does  in  tracers  (Jenkins  1982).  The  density  range  <r„ 
=  26.7-27.1  ( Fig.  4 )  is  about  1000  km  from  its  winter 
outcrop  (McDowell  et  al.  1982:  Keffer  1985).  Unlike 
internal  waves,  Ertel  vortieity  anomalies  cannot  prop¬ 
agate  because  they  are  tied  to  water  parcels.  They  can 
only  be  carried  into  the  ocean  interior  by  relatively 
slow  low-frequency  advection.  Because  of  weak  dissi¬ 
pation  in  the  interior,  anomalies  with  A_-  >  100  m  might 
survive  relatively  intact  as  they  are  subducted  and  swept 
around  the  gyre.  Elowever.  if  the  arguments  about  ero¬ 
sion  by  dissipative  turbulence  already  presented  are 
correct,  detrainment  cannot  explain  the  50-m  wave¬ 
length  anomaly  at  400-m  depth  as  the  measurements 
were  collected  six  months  after  retreat  of  the  winter 
mixed  layer,  twice  the  decay  time  quoted  above.  Given 
the  uncertainties  in  the  effects  of  turbulence,  a  more 
convincing  argument  might  be  that  buoyancy  forcing 
in  late  winter  would  create  0.  S  anomalies  as  well  and 
these  are  absent  between  depths  of  150  and  600  m. 
Therefore,  detrainment  in  late  winter  can  be  ruled  out. 

Generation  of  Ertel  vortieity  anomalies  on  the  flanks 
of  Ampere  Bank  or  other  nearby  seamounts  to  the 
northeast  (Fig.  I  )  followed  by  shedding  or  injection 


into  the  pycnocline  (Armi  1978;  Geyer  and  Signell 
1990:  Prater  1992)  is  the  most  plausible  formation 
mechanism  on  the  grounds  of  both  proximity  of  Am¬ 
pere  Seamount  and  the  intensity  of  irreversible  pro¬ 
cesses  at  topography.  GL^erved  dissipation  rates  in 
stratified  waters  on  the  Hanks  and  summits  of  sea¬ 
mounts  are  100-1000  times  greater  than  typical  pyc- 
nocline  dissipation  rates  (Osborn  1978:  Nabatov  and 
Ozmidov  1988;  Padman  and  Dillon  1991;  Toole  and 
Schmitt  1992.  personal  communication),  reducing 
formation  times  to  fractions  of  a  day. 

Evidence  for  benthic  boundary-layer  detachment 
near  topography  has  been  reported  by  Armi  (  1978). 
Osborn  (  1978 ).  Eiogg  et  al.  (  1978).  and  Nabatov  and 
Ozmidov  (  1988).  Rossbv  numbers  of  at  least  0.3  are 
revealed  by  the  drifter  tracks  (  Fig.  3).  Laboratory  ex¬ 
periments  in  a  stratified  fluid  produce  vortex  shedding 
for  Rossbv  numbers  greater  than  0.03  ( Zhang  and 
Boyer  1993).  Elowever.  the  numerical  simulations  of 
Smolai  kiewicz  and  Rotunno  ( 1 989 )  demonstrate  that 
inviscid  flow  can  produce  “vortices"  in  the  wake  of  a 
seamount  that  are.  in  fact,  manifestations  of  nonlinear 
internal  waves  with  no  Ertel  vortieity  signal,  raising 
doubts  about  interpretations  of  "eddy-shedding”  lab¬ 
oratory  experiments.  Bormans  and  Garrett  (1989) 
suggest  that  a  necessary  criterion  for  flow  separation  is 
that  L' /(Jr)  >  1.  where  r  is  the  radius  of  curvature  of 
the  topography.  For  the  observed  flow  speeds  of  5-10 
cm  s  l.  flow  separation  would  require  a  radius  of  cur¬ 
vature  less  than  1  km.  Features  of  this  size  can  be  seen 
on  the  200-  and  500-m  isobaths  in  Figs.  2  and  3. 

In  conclusion,  within  a  few  kilometers  of  Ampere 
Bank,  internal  waves  and  geostrophic  fluctuations  co¬ 
exist  on  wavelengths  A;  =  50-400  m  and  A/,  =  7-15 
km  (see  also  Kunze  1993).  Traditionally,  these  scales 
have  been  interpreted  as  internal  waves.  Vortex 
stretching  (strain)  and  available  potential  energy  are 
dominated  by  geostrophy.  while  both  internal  waves 
and  geostrophic  motions  contribute  to  the  horizontal 
kinetic  energy.  Despite  these  results,  it  is  premature  to 
assign  a  significant  contribution  from  Ertel  vortieity- 
carrying  dynamics  to  pycnocline  submesoscale  and 
finescale  fluctuations  because  the  most  plausible  source 
for  the  observed  Ertel  vortieity  variance  is  flow  inter¬ 
action  with  local  topography.  The  strength  of  Ertel 
vortieity-carry  ing  fluctuations  in  the  main  pycnocline 
away  from  topography  remains  to  be  determined.  In¬ 
direct  evidence  for  subinertial  shear  comes  from  hor¬ 
izontal  diffusion  in  the  ocean  being  much  higher  [  ~3 
nr  s  1  ( Ledwell  et  al.  1993)]  than  can  be  accounted 
for  by  internal  wave  shear  dispersion  [~0.01  nr  s  1 
(Young  et  al.  1982)].  Probing  the  dynamics  of  the 
submesoscale  and  finescale  promises  a  rewarding  area 
for  research  in  the  coming  decade. 
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S[^](A-)  =  E„f2  ~~  h2(x  -  arctan(X)) 

A  (A_-  +  A*)* 

( A4 ) 

/  x  \  A2A* 

x(arctan(x)-^T7)^TT^'  (A5) 

where  the  cutoff  aspect  ratio  X  =  (A  knjfk:)- 


APPENDIX  A 


The  Low-Pass  Filtered  GM  Spectrum 


The  GM  internal  wave  model  spectrum  (Garrett  and 
Munk  1979;  Munk  198 1 )  is  usually  posed  in  terms  of 
intrinsic  frequency  a >  and  vertical  wavenumber  k:  (or 
mode  number  j)  and  is  separable  in  these  variables. 
For  comparison  with  the  XCP  survey  data  in  Figs.  12 
and  13.  the  model  spectrum  needs  to  be  expressed  in 
terms  of  horizontal  and  vertical  wavenumbers,  then 
low-pass  filtered  in  horizontal  wavenumber  to  elimi¬ 
nate  the  high  horizontal  wavenumbers  that  are  sup¬ 
pressed  in  making  horizontal  least-squares  fits.  This  is 
particularly  important  for  the  horizontal  divergence  V 
=  (k/f  V„)and  relative  vorticity  f  =  (k„  X  V„),  which 
are  not  red  in  horizontal  wavenumber  but  have  much 
of  their  variance  at  small  horizontal  scales. 

Using  the  approximate  internal  wave  dispersion  re¬ 
lation  valid  for  w  <  A’,  which  should  apply  for  the 
horizontal  scales  resolved  by  the  survey. 


The  GM76  spectra  (Cairns  and  Williams  1976)  for 
horizontal  divergence  and  relative  vorticity  can  be  ex¬ 
pressed  as 


S[V](A„,  k:) 


EJN2N,J)2 


U2k ?  +  A '2A?,)  (k:  +  A?)2 


Sma„,  k:)  =  EJyN2N,Jh2 
k2n 


(j2k\  +  N~k2i)~  ( k:  +  A?)2 


where  the  dimensionless  energy  level  E„  =  6.3  X  I0“5, 
a  typical  upper-pycnocline  buoyancy  frequency  N„ 
=  5.2  X  10~3  s1,  b  =  1300  m,  and  the  scale  wave- 
number  A *  ~  ( x / b)(  N /  N„)  j *.  where  j *  =  3  at  mid¬ 
latitudes.  Integrating  ( A2 )  and  ( A3 )  from  A//  =  0  to  a 
high-wavenumber  cutoff  ku<  provides  a  low-pass-fil¬ 
tered  spectra  in  terms  of  vertical  wavenumber  A;  and 
the  cutoff  horizontal  wavenumber  A//r, 


APPENDIX  B 

Low-Pass  and  Contamination  Response  Functions 
for  a  Cross-shaped  Array 

An  array  containing  a  finite  number  of  elements 
cannot  determine  the  horizontal  gradient  exactly.  Such 
an  array  low  passes  scales  smaller  than  the  element 
separation,  but  imperfectly,  so  that  there  is  some  con¬ 
tamination  by  higher  wavenumber  variance.  In  partic¬ 
ular,  for  a  horizontal  array,  relative  vorticity  will  be 
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Fig.  Bl.  The  low-pass  attenuation  A' (solid)  and  contamination 
G  (dotted)  response  functions  for  a  cross-shaped  array  containing 
12  elements  per  leg  as  a  function  of  normalized  horizontal  wave- 
number  brk„,  where  hr  is  the  separation  between  adjacent  elements. 
The  low-pass  attenuation  response  function  is  one  for  hrk„  0.2 
and  falls  off  steeply  for  higher  wavenumbers.  The  contamination 
function  G  has  a  peak  of  0.03  at  hrk„  =  0.7. 
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Fit; .  B2.  Vertical  wavenumber  spectra  of GM  model  horizontal  divergence  ( solid )  and  relative 
vorticity  (dot)  filtered  with  the  low-pass  attenuation  and  contamination  response  functions  for  a 
cross-shaped  array  with  12  elements  per  leg  and  hr  -  0.2  km.  The  spectra  are  displayed  in  both 
log-log  ( left )  and  variance-preserving  ( right )  formats.  The  filtered  vorticity  spectrum  exceeds  the 
divergence  spectra  by  —  10'r-  for  =  20-50  m. 


contaminated  by  horizontal  divergence  and  vice  versa 
(Lien  and  Muller  1992b).  Lien  ( 1992.  personal  com¬ 
munication  )  has  formulated  the  low-pass  attenuation 
and  contamination  array  response  functions  for  a  cross- 
shaped  survey.  His  formulation  is  reproduced  here  with 
his  kind  permission.  The  low-pass  attenuation  response 
function  is 


F(k„)  = 


l _ 6. 

\M{M2  - 


)2  a/  at 

XH2J-M-1) 

j.  I  A- I 


X  (2k  -  M  -  1 )  x 
-  Mik  -j)k„br) 


M(k  -  j)k„br) 


-sin(2arcta„(- 


XJ2((  V( 2 7  -  M  -  \)2  +  (2k  -  M  -  \  )2)k„br/2) 

(Bl) 


and  the  contamination  array  response  function 
G{k,,)  ~  2{m(M2  -  dmJ 

M  M 

x  z  2  (2 j-M-  \  )(2k  —  M  —  1) 

y=l  A  -  I 

X  Jo((k  —  j)k/,br)  —  F(k/i),  (B2) 


where  M  is  the  number  of  elements  in  each  leg.  k,t  the 
horizontal  wavenumber  of  interest,  and  hr  the  sepa¬ 
ration  between  adjacent  elements. 

The  two  filter  functions  are  displayed  as  functions 
of  normalized  horizontal  wavenumber  k,fbr  in  Fig.  B! 
for  M  =  12.  The  low-pass  function  F  ( solid  curve)  is 
one  for  k„br  <  0.2.  and  falls  off  steeply  for  higher  wave- 
numbers.  The  contamination  function  G  rises  above 
the  low-pass  response  function  and  peaks  at  G’max 
=  0.03  for  slightly  higher  k,,br  =  0.7  so  that  relative 
vorticity  ( horizontal  divergence )  estimates  will  be  con¬ 
taminated  by  horizontal  divergence  ( relative  vorticity ) 
from  smaller  wavelengths. 

Figure  B2  displays  the  filtered,  contaminated  GM 
spectra  for  horizontal  divergence  and  relative  vorticity 

sm(k:)  =  J  \  SF](k:,k„)-F(k„) 

+  Sin(k;,k„)-G(kfl)\dk„ 

S[V](k:)  =  J  {SU](k:.  k„)- F(k„) 

+  S[V](k:.kll)-G(k„)}dkll  (B3) 

with  M  =  14  and  br  =  0.3  km  appropriate  for  survey 
1 .  The  vorticity  spectra  will  be  more  contaminated  by 
divergence  than  the  divergence  by  the  vorticity  because 
the  GM  divergence  spectra  contain  more  variance  at 
high  ho.izontal  wavenumber.  The  filtered,  contami¬ 
nated  GM  spectra  differ  from  the  simple  boxcar  filter 
used  for  Fig.  1 3  only  in  that  the  vorticity  exceeds  the 
divergence  spectra  by  10%  at  high  vertical  wavenum- 
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bers  ( X.  <  60  m).  We  caution  that  a  cross-shaped  survey 
is  a  suboptimal  sampling  pattern,  equivalent  to  four 
points  on  a  circle  (Lien  1992.  personal  communica¬ 
tion).  and  succeeds  in  this  paper  only  because  the  hor¬ 
izontal  divergence  and  relative  vorticity  are  of  com¬ 
parable  magnitude.  For  smaller  surveys  sampling  the 
scales  for  which  internal  wave  horizontal  divergence 
greatly  exceeds  the  relative  vorticity.  the  optima!  array 
pattern  would  be  a  circle  or  a  square  grid  containing 
at  least  nine  stations. 
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